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Deformation of Ising Hamiltonian by means of replacing a site spin s, by s* and statistics gen- 
eralization with help of the substituting deformed probability p\ instead of pi are studied jointly 
within mean-field scheme. Such deformed model is shown to be related to the phase transition 
of the second order with unusual set of critical indices depending essentially on the deformation 
parameter q. Scaling relations turn out to be invariant with respect to the deformation. 
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I. INTRODUCTION 



II. MAIN STATEMENTS 
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The Ising model is known to be a corner stone of the 
microscopic theory of phase transitions Q — • Evident 
simplicity of this model is based on a supposition that ev- 
ery lattice site i has a spin taking only two magnitudes 
S{ = ±1. Although exact solution is not found for the 
three-dimension Ising model, its using allows one to ex- 
plain main peculiarities of real phase transitions. Along 
this line, the qualitative picture becomes clear already 
within framework of the mean-field approach. 



Recent considerations pay much attention to study 
complex systems which self-similarity derives to general- 
ization of the Gibbs-Boltzmann statistics to Tsallis-type 
one (see Ref.^jJ and references therein). Formally, such 
a generalization is performed by means of replacement 
of the probability p t by the so called deformed probabil- 
ity pj with a positive index q < 1. Here, we propose to 
complete such a procedure by the relevant deformation 
on the microscopic level. Namely, we deform the Ising 
Hamiltonian by means of replacing a site spin Si by sj. 
Such a deformation is shown to uphold the second or- 
der phase transition, whereas the set of critical indices 
becomes depending on the parameter q essentially. 



The Letter is organized in the following manner. In 
Section 2 the q deformed Hamiltonian is postulated and 
simplified within the mean-field approach. Section 3 de- 
als with determination of the fractional average allowing 
us to obtain the definition of the order parameter. In 
Section 4 q-deformed distribution function is found to 
write a self-consistency equation for the order parameter. 
Section 5 devotes to the consideration of asymptotic so- 
lutions of the self consistency equation. The last Section 
6 contains a brief discussion concerning scaling relations. 



We postulate the q-deformed Hamiltonian in the fol- 
lowing form: 



(1) 



Here, summation runs over N lattice sites i ^ j, Jij is 
effective interaction potential, h is external field, Sj = ±1 
is a site spin, q < 1 is a deformation parameter. 

In the simplest way, we use further the usual scheme of 
the mean-field approximation 0-0. Within this ap- 
proach, one follows to replace the multiplier by an av- 
erage value (s q ). Moreover, we shall take into account the 
only interaction of nearest neighbor sites whose number 
equals z and potential is reduced to the constant J > 0. 
Then, the mean-field effective Hamiltonian is as follows: 



H, 



ef 



h q s1, h q = h + T c (si) (2) 



where a characteristic temperature T c = zj is intro- 
duced. 

It is easy to foresee the postulated Hamiltonian will be 
meaningless if the index q < 1 is not chosen to satisfy 
the condition 



1. 



(3) 



Then, the complex representation — 1 = c J7r ( 2 ™+ 1 ) ; n 
0, ±1, . . . arrives at the set of rational numbers 



<l 



2m + 1 
2n+ 1 : 



m < n 



(4) 



with integers m,n = 0, ±1, Obviously, above reduc- 
tion of the continuous set q £ [0, 1] into the manifold of 
rational numbers 10} is caused by the discrete symmetry 
of the Ising model being invariant with respect to the 
transformations {si} — > {— Sj}, h — > —h. 

As a result, the site energy in the effective Hamiltonian 
(j2J takes the linear form 
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III. CALCULATION OF THE FRACTIONAL 
AVERAGE 

Now, we should express the fractional average (s q ) by 
means of the order parameter r\ = (s) . Obviously, such a 
problem will be quite correct only for self-similar systems 
where it is reduced to definition of an index p(g) in the 
relation p| 



(s q ) = rf {<l \ 



(6) 



If relevant distribution varies very slightly near the ori- 
gin x — and decays abruptly in the limit x — > 00 (see 
Eas. l|ll[) . (|13|l below), we have the following estimations: 



These derive immediately to the principle relation 



(X*)~{X)**\ p(q) 



(8) 



Respectively, the effective field in Eq.JSJ) takes the form 



h q =h + T c r)— 



(9) 



parameter 77 by the power-law function 77 



usual manner the set of g-deformed hyperbolic functions: 
sinhq(a;), cosh 9 (a;), tanh g (ar) and coth 9 (x) Q. 
We calculate now the site partition function 

z q = exp g (-/3<l) g£i )=2cosh ? (/?<l)A) (14) 

Si=±l 

that differs from the usual one with g-deformation and 
(l) q factor appearance. According to Eas. i(TT)l . 
the latter is determined by the equation 



(1),= 



[exp g {(3{l) q h q )] q + [ex Pg (-(3(l) q h q )] q 
[ex Pg (/3{l) q h q ) + exp q {-(3{l) q h q )} 9 



(15) 



With accounting Eo.l|14jl. the order parameter rj = (s) is 
defined through the self-consistency equation 



i+a 

T] 2 = tanh 9 



(i>Afa)] 



(16) 



that differs from the usual one with change of the order main 



following from Eqs.©, lfTT|l . (fT^I and condition J2J). 

Solutions of the system l)15[l. (|16|l are shown in Fig.l 
for the external field h = and different values of index q. 
Main peculiarity of the temperature dependencies T](T), 
(l)q(T) is in decreasing the order parameter < r\ < 1 
accompanied by the relevant increasing the parameter 
1 < < 2 at the temperature growth within the do- 



T cq <T<T q , T cq = (l-q)T c , T q = 2 1 - q T c (17) 



IV. STATISTICAL SCHEME 

Along the line of using g-deformed quantities, we need 
to state on the Renyi entropy 



^(l-^-Mn^P^} 



(10) 



which is easy seen to take the Boltzmann form in the 
limit q — > 1. Then, the maximum entropy princi- 
ple taken with accounting the normalization condition 
53/ Si } P{si} = 1 and the definition of the g-deformed 
internal energy E q = J2{ Si } £ i{ s i)P q { s i} derives to the 
distribution function 



{i) q = E pKsi) 



(12) 



where z q is the specific partition function related to one 
site, /3 = i is inverse temperature measured in the energy 
units; g-deformed exponential is as follows: 

expJx) ee { t 1 + (! - at 1 + (1 - q)x > 0, 

q [0 otherwise. 

(13) 

Such a form is known to generalize usual exponential 
related to the limit q — > 1. It allows one to generate in 



(at T < T cq one has 77 = (l) q = 1, whereas at T > 
T q there are r\ = 0, (l) g = 2). Thus, decreasing the 
deformation parameter 1 > q > derives to monotonous 
growth of the order parameter that transforms smoothly 
falling down dependence r](T) into step-like one. 



Fig.l. Temperature dependencies of both the renor- 
malization parameter (upper panel) and the order pa- 
rameter (low panel) for different indexes q (curves 1 — 7 
relate to q = 0, 1/17, 1/7, 3/11, 5/11, 5/7, 1, respec- 
tively). 



As demonstrated in Fig. 2, the effect of the external 
field is similar to the usual case q = 1: high-temperature 
dependence rj(h) has the monotonically growing form 
(dashed curves) that takes, at the low temperatures 
T < T q , the falling down domain located in the vicinity 
of the point h = 77 = (solid curves). At fixed tem- 
peratures, decrease of the deformation parameter is seen 
to arrive at the more sharply defined non-monotone de- 
pendence. At critical field h q = — T c (l — T/T cq ), the 
dependence rj(h) undergoes a break, whereas in the limit 
of the small fields h <C h q one has 



h/T c 



1 - T/Tq 



(18) 
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Fig. 2. Field dependencies of the order param- 
eter: dashed curves 1 6 correspond to the 
fixed index q — 5/11 and temperatures T/T c = 
1.5, 1.7, 1.9, 2.2, 2.5, 2.8; solid curves 1 — 6 relate to 
indexes q = 0, 1/7, 3/11, 5/11, 11/17, 1 at T = 0.5T C . 
Insertion shows details of the coordinate origin neighbor- 
hood. 



V. ASYMPTOTIC REGIMES 

To obtain analytical description, we look over the pos- 
sible asymptotic solutions of the system ()15|l , (|16|) in lim- 
iting cases of both zero and non-zero external fields. 



critical temperature T q . Here, the hyperbolic cosine in 
Eq.lfHfl can be expanded over the argument (3(1) q h q ~ 

pT q r) 2 in series comprising of even terms only. In accor- 
dance with relevant temperature dependence (|21(l . these 
terms have orders [(T q — T)/T q ] n with the lowest index 
ri = 2 because the first-order term (n = 1) is suppressed 
by self-action effects. Thus, in the limit (T q -T)/T q < 1 
we obtain F ~ T q [(T q - T)/T q ) 2 and the specific heat 
C ~ d 2 F/dT 2 does not vary with temperature. This 
results in the magnitude 

a = (27) 

of the critical index of the dependence C ~ \(T q — 
T)/T q ]-« 



Zero field 



Non— zero field 



For small order parameters (77 <C 1), one has the ex- 
pansions 



tanhq (x) 



g(2-g)„3 



(1) 



2 l-q 



l + 2 2(i-«)(i_ g )2 ^ hq y 



(19) 
(20) 



Then, the temperature dependence of the order parame- 
ter takes the power-law form 



A : 



(3-2q)(2 9 -l)' 







1+9' 1 i 



2 1 ~iT r 



(21) 
(22) 



that corresponds to critical domains T q — T <C T q in 
Fig.l. At q < 1, the index (3 exceeds the usual magnitude 
0=1/2. 

At the marginal magnitude q — 0, one has 



tanho(x) = x, (l)o = 2 
and equation (|16fl arrives at the condition 

T 



(23) 



(24) 



It remains valid for arbitrary values 77 at T — T q . So, this 
case corresponds to the step-like curve 1 in Fig.l. 
In the limit q <C 1, there are the expansions 



(l) q ~2 1 - q , tanh g (x) ~ x ( 1 - 2q 



1 - x< 



Respectively, the order parameter 

2q 



V 



1 



1 - T/T q 



(25) 
(26) 

q «: 1 



decays very fast within tight window 1 — T/T q 
(see curve 2 in Fig.l). 

Finally, we consider the form of the temperature de- 
pendence of the free energy F = —TN\n(z q ) near the 



Here, we start with the approximate equation 

2 2(l-g) 



2 1 - q f3h q 



1 



-iflKY 



(28) 



following from Eas. (|16|l . (|19|l (notations h q and A are 
given with Eqs.®, J22J)- In the limit h — ► 0, differen- 
tiation of Ea. ((28)) with respect to the field h yields the 
following expression for the susceptibility \ = drj/dh: 



T cX = 2(l + q)- 1 V - 



T-T q 



■,1+9 



A 



(29) 



Thus, in disordered state (77 = 0, T > T q ) one has the 
susceptibility % = 0. Physically, this means a suppression 
of critical fluctuations in the disordered phase — quite 
differently from the usual picture. 

In ordered phase (rj ^ 0, T < T q ), inserting Eq.JSTJ 
into Eq. (|29|l we get the power-law dependence 



,1H 

A 2 r+3 



T q -T 



with the critical index 



7 



1 1 

l -2l 



q 



(30) 



(31) 



With decreasing the deformation parameter 1 > q > 0, 
this index falls down monotonically from the usual value 
7 = 1 to 7 = 1/2. 

In the opposite limit h — > 00, Ea. (|28fl related to the 
critical point T — T q derives to the power-law depen- 
dence 



with the critical index 



5(1 + *)- 



(32) 



(33) 
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VI. DISCUSSION 

Above mean-field consideration shows the deformation 
of Ising Hamiltonian by means of replacing a site spin 
by a power function sf derives to the phase transition 
of the second order with unusual set of critical indices 
(|77jl . [|3T1) and gg. The zeroth magnitude of the 

first of them is caused apparently by the mean field ap- 
proach. It is easy to convince also in zeroing indexes 
e, £ defined by the field dependence of the specific heat 
C ~ hr e and the space dependence of the correlation 
function G(r) ~ r~( d ~ 2+l ») Moreover, there are crit- 
ical indexes fi and v defining the field-temperature de- 
pendencies £ ~ and £ ~ (T q — T)~ u of the correlation 
length £. Among complete set of these indexes, only two 
are independent, whereas the rest are known to be de- 
termined by the scaling relations. Then, the question of 
fundamental importance arises how influences the defor- 
mation on these relations. 

The first of them, Widom relation, takes the usual form 

□ 

(36 = (3 + 7 (34) 

that follows from comparison of two chains of definitions: 
rj ~ x h^ h ~ (T q - TY+^i and h ~ rf ~ (T q - Tf s . 
Using the relation £ ~ /i _A1 in the latter, one has 

H(fi + l) = v. (35) 



Combination of Eas. (|34() . I|35(l yields the relation (36 fj, = 
v that together with Eas. (|22l) and (|33|) gets 3^ = 2v. 
Finally, using the fluctuation-dissipation theorem derives 
to the equality @ 

2v = 7 (36) 
that defines the rest of indexes: 

What about the Essam-Fisher equality, it takes the 
form 

a + (3 + 7 = ~ (38) 

following from the consequence C(T q — T) 2 ~ rj~rh =>■ 
h ~ (T q — T)i~ a that differs from the usual one by re- 
placing rj by fjT 1 . 

It is worthwhile to note all of equalities l|34l) - 13611 . 
(|38|l are invariant with respect to the deformation. 
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